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The ARCADE 2 collaboration and other experiments have reported a significant excess in the 
isotropic radio background, whose homogeneity cannot be reconciled with clustered sources. This 
suggests a cosmological origin prior to structure formation. We investigate several potential mech- 
anisms and show that injection of relativistic electrons through late decays of a metastable particle 
can give rise to the observed excess radio spectrum through synchrotron emission. However, con- 
straints from the cosmic microwave background (CMB) anisotropy, on injection of charged particles 
and on the primordial magnetic field, present a challenge. The simplest scenario is with a > 9 GeV 
particle decaying into e + e~ at a redshift of z ~ 5, in a magnetic field of ~ 5/iG, which exceeds 
the CMB _B-field constraints, unless the field was generated after decoupling. Decays into exotic 
millicharged particles can alleviate this tension, if they emit synchroton radiation in conjunction 
with a sufficiently large background magnetic field of a dark U(l)' gauge field. 



1. INTRODUCTION 

The ARCADE 2 collaboration has measured the dif- 
fuse radio background at several frequencies between 3 
and 90 GHz, with relatively small errors up to 10 GHz [1]. 
Their measurements extend earlier ones at lower frequen- 
cies 0.01 — 1 GHz [2]- [5] that arc well fit by a power-law 
dependence of the brightness temperature on frequency, 

T ~ T R {v/v Q f (1) 

with (3 = — 2. 6± 0.036. The measured value of j3 is within 
the ballpark of expectations (/? = —2.7) for sources in 
which synchrotron emission dominates (e.g., star-forming 
galaxies), but the amplitude Tr = 1.26 ± 0.09 K (with 
v o = 1 GHz) is ~ 6 times too large to be explained by 
extrapolations of populations of known resolved sources 
[6] . Attempts to identify the excess with previously over- 
looked standard astrophysical sources have so far not 
met with success [7]- [9]. There have been several at- 
tempts to provide a new origin through the annihilation 
of dark matter into charged particles, which subsequently 
undergo synchrotron emission [10]- [14]. It was recently 
pointed out that there is a serious challenge to such mech- 
anisms: the spatial fluctuations in the observed radio sig- 
nal are much too small to be consistent with sources that 
have undergone clustering [15], as would be the case for 
dark matter annihilating near the present time. 1 

These considerations suggest a cosmological origin of 
the radio excess, from an epoch prior to the formation 
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1 Ref. [14] studied the possibility of dark matter annihilating into 
charged particles at earlier times, assuming that ultracompact 
minihalos and very large magnetic field exist, finding that Comp- 
ton scattering leads to overproduction of the diffuse x-ray back- 
ground relative to observed values. 



of structure and consequent large inhomogeneities in the 
density of decaying dark matter or of the cosmological 
magnetic fields. On the other hand, the excess photons 
should have been produced relatively late, after redshift 
z = 1100 when recombination occurred, to avoid ther- 
malizing away such a spectral feature. When one tries to 
imagine a mechanism for producing excess photons with 
the right spectrum at this epoch or later, it proves to be 
highly constrained, due to cosmic microwave background 
(CMB) constraints B < 10~ 9 G on primordial magnetic 
fields and on injection of ionizing charged particles. The 
purpose of the present work is to point out some mech- 
anisms that could be promising and the challenges that 
they must be overcome. We show that a spectrum of 
excess diffuse radio background consistent with observa- 
tions can be generated either through synchrotron emis- 
sion from electrons injected into the plasma at some time 
after decoupling, depending upon the mass of the decay- 
ing particle that produces them. However, CMB con- 
straints on the magnetic field (if it was generated before 
decoupling) and on the amount of injected energy turn 
out to be in conflict with this hypothesis. 

For example, we find that a metastable particle x that 
decays via x e + e~ can produce the desired radio back- 
ground through synchrotron emission if m x > 700 McV, 
consistent with the CMB constraint on ionizing radia- 
tion, but it requires a magnetic field some 100 times 
greater than the CMB constraint on primordial B fields, 
necessitating a late-time magnetogenesis mechanism. A 
more complicated model may be viable: if x decays into 
millicharged particles that also carry a dark U(l) gauge 
charge, and if there is a corresponding background of 
dark photons with a sufficiently large dark magnetic field, 
the CMB constraint on charged particle injection can be 
robustly overcome. 

Our paper starts with a brief recapitulation of the ob- 
served radio spectrum, reviewing in section 2 the data 
and the possibility that a spectral index —2.5 for the ex- 
cess temperature provides a consistent description. In 
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section 3 we note that direct production of photons 
through decays or annihilations does not give this kind 
of spectrum, nor do processes involving low-energy elec- 
trons, motivating the alternative of relativistic charged 
particle injection. Section 4 discusses the rate of en- 
ergy loss of charged particles through Compton scatter- 
ing (CS) and synchrotron emission, and in it we derive an 
upper bound on the primordial magnetic field around the 
decoupling epoch, from CMB constraints. This bound 
implies that CS is by far the dominant means of en- 
ergy loss. Section 5 computes the spectrum of Compton- 
scattered photons and consequent bounds from x-ray ob- 
servations. Section 6 shows that synchrotron emission 
can produce a radio background with the desired proper- 
ties, if the initial electron energy is sufficiently large. In 
sections 5-7 we make a simplifying assumption that the 
charged particles are injected at specific time (sudden 
decays). In section 7 we show that this approximation 
works very well for decays, and we generalize it to the 
case of annihilations. In section 8 we show that CMB 
constraints arc in conflict with the above mechanisms if 
the injected particles are electrons, and if the primordial 
B field was generated before decoupling. In section 9 
we suggest a more exotic scenario in which millichargcd 
particles are instead injected, as a second example of how 
the challenges might be overcome. Conclusions are given 
in section 10. A possibly novel way of generating primor- 
dial magnetic fields (though too weak for our purposes) 
is elaborated in appendix A. In appendix B we show 
that brehmsstrahlung cannot yield the observed spectral 
shape. 
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FIG. 1: Excess antenna temperature T A versus frequency as 
measured by refs. [1]- [5]. Dotted and dashed lines show u~ 2S 
and u~ 2 ' 5 power laws, respectively. 
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2. THE OBSERVED SPECTRUM 

Taking data from table 4 of ref. [1], which also com- 
piles results from the earlier experiments [2]- [5], we have 
rcplottcd the antenna temperature 2 T A data versus fre- 
quency and display them in fig. 1. The ARCADE-2 data 
are the 6 highest frequency points (we do not show those 
above 10 GHz, where the error bars are much larger). 
Also shown are power law dependences, with the 
ARCADE-2 best fit value (3 = -2.6, and the theoreti- 
cal prediction ft = 2.5 that we will derive in section 6. 
From fig. 1, it is clear that the preference for ft = —2.6 
is driven by the measurements near v — 45 MHz and 3.3 
GHz that have the smallest errors. This assumes there 
is no systematic miscalibration between the two different 
experiments [3] and [1] that determined the respective T A 
values. Even if not, -2.5 is within 2.8<r of ARCADE-2's 
best-fit value. Prior to the ARCADE-2 measurement, 
j3 = —2.5 gave a good fit to the existing data, as was 



2 The antenna temperature takes the place of E/{e@ E — 1) in the 
familiar formulas for fluxes, intensities, etc., of blackbody radia- 
tion. We use units h = c = kg = 1. 



FIG. 2: Temperature of cosmic background radiation to 
higher frequencies/energies, compared to E~ 2 ' 6 and E~ 25 
power laws, data taken from fig. 1 of ref. [16]; lines are our 
fits to the current low-frequency observations. Region of fig. 
1 is shown by box in upper left-hand corner. 

pointed on in ref. [2]. In the present work, we will con- 
sider /3 = —2.5 to be adequate for describing the obser- 
vations. It will be seen that there are more challenging 
problems for finding a working mechanism than the dis- 
crepancy between —2.5 and —2.6. 

To put the radio observations into a broader perspec- 
tive we have taken the compendium of diffuse radiation 
backgrounds presented in fig. 1 of [16] and translated it 
into antenna temperature using T A = 47r 2 ^ _2 /„ to in- 
dicate how a new source in the radio region (energies 
less than ~ 10 -6 cV) would compare to higher-energy 
backgrounds if it continued with the same power law in- 
definitely. At high energies E > 10 cV, the index —2.5 
or —2.6 would not be sufficiently steep to remain within 
observational bounds, but we will see that this need not 
be a serious limitation. First, the mechanism we pro- 
pose has a natural cutoff energy above which the spec- 
trum of the new contribution drops sharply, and this can 



be dialed by adjusting the initial energy of the injected 
charged particles. Even in the absence of such a cutoff, 
one must consider whether the universe is transparent 
to radiation in a given energy range and at a given red- 
shift. The transparency window as a function of z and 
Ej is given in refs. [17, 18]. At low rcdshifts, this window 
extends from roughly 1 keV to 10 5 GcV, and it is also 
open at E < 10 eV below the energy of Lyman- a ab- 
sorption. Thus any x-rays produced in conjunction with 
the radio excess will be absorbed by the baryonic plasma 
for energies below ~ 1 keV. At higher redshifts this cut- 
off increases to ~ 10 keV. Depending upon the details, 
this effect could weaken the constraining power of exper- 
iments like Chandra [19] and XMM-Newton [20] that we 
will discuss in section 5. 



3. MOTIVATION FOR RELATIVISTIC 
CHARGED PARTICLE INJECTION 
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FIG. 3: Comparison of the shape of the observed radio ex- 
cess flux (in arbitrary units) to that produced by dark matter 
decays \ ~ ► 



One might hope to bypass the complications of syn- 
chrotron and Compton emission by a mechanism that 
directly produces photons after decoupling. An axion- 
like particle \ with a very long lifetime and which decays 
into two photons is the simplest possibility. Although 
it produces monochromatic photons at any given time, 
their energies are rcdshiftcd and so a continuous spec- 
trum results. Rcf. [21] finds, in the case where the de- 
caying particle lifetime r x is greater than the age of the 
universe, that the flux of such photons is given by 
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where A = 10- 7 (cm 2 s str)" 1 ^ 1 ^ /r x ) (10 GeV/m x ) 
assuming \ is the dark matter, and k = U,^/U, m = 3. 
The energy-dependence can be derived from the integral 
over redshifts [22] up to that at decoupling (z^ oc ) 
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where the Hubble parameter is proportional to 
\/(l + z) 3 + hi and the delta function yields the spectrum 
of photons from a given decay. 

On the other hand, the radio signal of interest has the 
observed flux 
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where T(E 7 ) ^ T R ■ (E y /E Q )~ 2 - 6 with E Q = A A x 10~ 
eV (corresponding to v = 1 GHz) and Tr = 1.1 x 10~ 



eV from (1). The shapes of (2) and (4) are compared in 
fig. 3. The slope of the decay spectrum is never as neg- 
ative as the observed one. We find that the agreement 
is not improved by considering decays into three pho- 
tons, approximating the spectral distribution as a box in 
place of the delta function in (3), nor by decreasing the 
lifetime such that the x' s disappear before the present 
time (which also changes the shape of the present-day 

spectrum via an extra factor e~^ to ^ Tx ^ 1+z ^ 3/2 in (3)). 4 
Neither do annihilations XX 77 produce better results. 
We conclude that particle physics mechanisms for direct 
production of photons do not work. On the other hand, 
the radio excess has a spectrum that appears to be con- 
sistent with synchrotron emission. Production of charged 
particles thus seems like a more promising approach. 

In order to generate the correct spectrum, we will fur- 
thermore require the injected charged particles to be rel- 
ativistic: E e m e . The cyclotron output from non- 
relativistic electrons in the cosmic magnetic field will be 
in the form of peaks around the gyrofrequency. This 
is at most ~ 1 Hz for values of B ~ nG allowed by 
CMB constraints — far below the MHz-GHz region of the 
ARCADE-2 excess. 

Compton scattering with the CMB is a second mech- 
anism which may produce low-energy photons. By com- 
puting the lower-limit to the scattered photon energy, we 
can show that this mechanism is not an efficient means 
of producing radio photons, if the injected electron is 
nonrclativistic. A CMB photon with energy E n i <C m e 
scattering head-on with a low energy electron of velocity 
n e < 1 follows the Compton formula, producing a photon 



3 Ref. [22] derives the spectral intensity /(A) rather than the flux, 4 To avoid the problem of large angular fluctuations in the radio 
giving an extra factor of A 2 /E~, ~ (1 + -z) 3 in the denominator excess, we should insist on lifetimes shorter than the time scale 

of the integrand. for structure formation 
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of energy E 1 : 
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(5) 



We are computing an upper limit to energy loss which we 
therefore maximise by choosing cos# = 0. This yields: 
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To explain the ARCADE excess, we are seeking a mech- 
anism to produce photons < 0.01-Ey (the average 
CMB photon energy). This is clearly not possible in the 
v e <C 1 regime of (6). 



4. INJECTION OF ELECTRONS AND 
MAGNETIC FIELD CONSTRAINTS 

We have argued that a more likely mechanism for pro- 
ducing the radio excess would involve injecting charged 
particles with excess energy into the primordial plasma 
at some time after recombination but before structure 
formation, and relying upon the subsequent synchrotron 
emission. For example, a metastable particle \ °f mass 
m x > 2m e could decay into e + e~ pairs at the appro- 
priate time. Alternatively a first order phase transition 
of some scalar field coupled to the Higgs boson could 
cause a small change in the Higgs vacuum expectation 
value, leading to a change in the mass of the electron 
and hence its kinetic energy. However it seems difficult 
to obtain relativistic electrons as required through the 
latter mechanism. The decaying x ~> e + e~ scenario on 
the other hand naturally provides relativistic electrons. 

It is enlightening to estimate how much energy is 
needed to produce the observed excess. Supposing that 
the spectral anomaly continues up to some maximum en- 
ergy .E max , its energy density is given by 
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(Recall that Eq corresponds to the frequency Uq = 1 
GHz.) As a fraction of the critical density, this is 

fir = 1.3 x HT 10 (£ max /10£ ) a4 . 

Eq. (7) gives a lower bound on the excess kinetic energy 
in electrons that needs to be injected. But in fact much 
more is required if we rely upon synchrotron emission, 
since it is not the most efficient means of dispersing the 
excess energy in the early universe; Compton scattering 
on CMB photons (7') is much more important. The total 
energy loss rate is due to the sum of the two processes, 
and their relative importance depends upon the energy 
density in CMB photons versus that in magnetic fields: 
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where 7 = E e /m e and a T = (87r/3)e 4 /m 2 is the Thom- 
son cross section. The relative efficiency for producing 
synchrotron radiation is therefore 
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Here we have alluded to the CMB constraint on the co- 
moving magnetic field B < O(10~ 9 )G [23]- [28]. Taken 
at face value, a limit of B < 1 nG combined with (7) 
implies that 10 7 times more energy in Compton photons 
is produced compared to synchrotron. If we wish to use 
the latter to produce the radio excess, it would then lead 
to an extra contribution to in photons via CS of order 
10~ 3 , which is 25 times greater than that in the CMB. 

To understand whether synchrotron emission can be 
more efficient than indicated by the above discussion re- 
quires us to examine the CMB contraint on B in greater 
detail. The primordial B field is believed to be stochastic, 
with a spectrum that in general depends on the length 
scale. The limit O(10 -9 ) G is thus a bound on the av- 
erage value B\ = (S(a;)i3(y)) 1 / 2 at some comoving scale 
A = \x — y\, which is usually taken to be 1 Mpc for 
CMB constraints. If the spectrum of B is scale invariant 
(meaning that n B = —3 in B\ ~ A~( ns+3 )), then there is 
no possibility that B\ could be larger than the nominal 
bound at some smaller scale. But if it has a blue-tilted 
spectral index n B > —3, then B\ could conceivably ex- 
ceed O(10 -9 ) G by some factor at sufficiently small scales 
(yet large enough to be coherent over the Larmor radius 
of an electron and thus produce synchrotron radiation) . 

A limitation is that magnetic fields are damped at 
small scales by radiative viscosity, and the damping scale 
is a function of the B field itself. The result of [24] for 
the damping scale can be expressed as [25] 
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(where n = n B ). This applies to vector perturbations, 
which were shown in rcf. [28] to dominate over scalar 
perturbations. At BiMpc = 10 _9 G, A D ranges from 0.05 
to 0.4 Mpc as n B goes from —3 to 2 (the maximum value 
usually considered for a causal mechanism of primordial 
magnetogenesis) . For a given n B , the maximum field 
value will be attained at the minimum undamped scale, 
just above X^. 

CMB constraints on B\ and n B are correlated, as 
shown in rcf. [27-29]. We reproduce the 95% C.L. limit 
on n B versus B\ from ref. [28] in fig. 4(a). The upper 
limit on B\ is strongest for the largest values of n B . We 
have evaluated B\ at the damping scale X D for points 
along this contour to find the maximum allowed B\ as a 
function of n B . The result, plotted in fig. 4(b), is that 
B\ D as large as 50 nG is allowed for a spectral index in 
the range n B = ( — 1.5, —1), corresponding to damping 
scales X D = (0.06 — 0.1) Mpc. Using this field strength 
in (9), the relative efficiency to produce synchrotron ra- 
diation is enhanced by a factor of 2500, and the energy 




density produced by associated Compton scattering be- 
comes around 0.01 of that in the CMB, if the excess radio 
signal is due to synchrotron emission. This relaxation of 
the nominal bound of B < a few nG is only useful if 
synchrotron emission occurs close to the decoupling era 
since thereafter we expect the damping scale to increase 
again as the free charge density decreases. This caveat 
will not affect our conclusions in the end since we will 
have to invoke either late-time magnetogenesis or other 
exotic particle physics to construct a working scenario. 

Stronger limits on B\ for scale-noninvariant spectra 
can be obtained by requiring that the energy density p B 
in B not exceed the critical density at earlier times [30]; 
the fraction of energy in ps increases at earlier times be- 
cause Xd depends upon time (eq. (10) is valid around 
the epoch of matter-radiation equality). However these 
limits would be evaded if either the time of magnetoge- 
nesis was not too early, or if there is an intrinsic cutoff 
in the initial spectrum (as one would expect there must 
be) that falls below the damping scale before p B becomes 
too large. Thus the upper limit -B max — 50 nG is fairly 
conservative, from the perspective of making the fewest 
assumptions about the unknown origin of the primordial 
magnetic field. Of course, if magnetogenesis occurs even 
later than recombination, we can evade this constraint as 
well. 

A single mechanism that could accomplish both of 
these at once might interesting for explaining the radio 
excess. For example, the decay process x ~~ ^ e + e~ pro- 
duces a small electric current at each decay, which gives 
rise to stochastic B fields. We can estimate their strength 
as 



(B(x)B(y)) 
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for relativistic electrons where n e is the electron density. 
In appendix A we derive this, and show that it is too 
small for our purposes. In the following, we will assume 
that magnetogenesis occurred before decoupling, so the 
maximum B field was no more than 50 nG, although the 
possibility of relaxing this assumption should be kept in 
mind if one is looking for loopholes. 



5. COMPTON SCATTERING 

Since we have argued that much more energy will be 
produced by Compton scattering than by synchrotron 
emission, it is necessary to compute the spectrum of CS 
photons resulting from injected electrons. Although this 
spectrum itself does not have the right properties to ex- 
plain the excess radio background, it has the potential 
to make a larger contribution to the diffuse x-ray back- 
ground than is observed, and thus provides a constraint 
on mechanisms that rely upon synchrotron emission to 
account for the radio excess. We derive such a constraint 
in this section. 

The rate of Compton emission of photons of energy E 1 
per scattered electron of initial energy E e 3> m e , involves 
an integral over the energy spectrum dn/dEy of CMB 
photons: 

f d E^^^f c ( q ,T)g c (E,/E Y ) 

(12) 
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where (following ref. [31]) 

f c (q,T) = 2glng + (1 + 2g)(l - g) + (13) 
with T = 4EyE e /m 2 e and q = E 1 /{T{E e - E 7 ) and 
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E 7 J Ej 
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The above parametrization is approximate, and can be 
inferred from fig. 3 of [31]; its main purpose is to show 
that the spectral index of Compton radiation changes by 
1 between the upscattering and downscattering regimes. 
/ c (g,r), which is meant to describe the regime where 
Ey > Ey , becomes roughly constant at its maximum 
value for E-y < Ey . This is true provided that the elec- 
trons are relativistic. The factor g c {E 1 / 'Ey) is designed 
to give the correct behavior for E 7 < Ey . We will further 
approximate (12) by replacing Ey with its average value 
Ey = 2.7 T<i, and dn/dEy = n 1 ^jEy where T<j is the 
temperature at time when the electron was injected (the 



6 



0- 



factor of z 2 d (by dimensionality of dn^/dE-y). Thus 



> 



■10- 



£ -20 



-30- 




2001 observations 
E ' power law 

E 2 6 power law 
Compton spectrum 



-4 -2 

log 10 E y (eV) 



FIG. 5: Solid (blue): the maximum allowed Compton spec- 
trum consistent with Chandra constraint, shown for m x — 
600 MeV. Shaded box shows lowest-energy Chandra observa- 
tion [19], near edge of transparency window below which x- 
rays are absorbed by intergalactic medium. Other data taken 
from fig. 1 of ref. [16] 



decay time td of x m the model where \ e + e~) and 
n 7:( 2 is the density of CMB photons at this temperature. 
Later on we will relax this "sudden decay approximation" 
by averaging over td with the appropriate dependence of 
n x on td- The spectrum of emitted photons at the decay 



time is then given by 
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(15) 

Wc ignore diffusion since the universe is homogeneous 
and isotropic during the epoch of interest, and we assume 
that the electron spectrum remains monoenergetic. The 
electron energy diminishes very quickly with time; by 
solving eq. (8) with p 7 > = p 7t d ^ Pb, we find 
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for t > td, where Ed = m x /2 is the initial electron energy 
when it was injected at time td, and 
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The time scale for energy loss is much smaller than the 
Hubble time, so we can ignore redshifting of the electron 
energy. It is convenient to change variable t — > E e = 
E e /E in the integral (15), using dE e /dt = —C e E 2 /E d , 
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The approximation 0.088 for the integral holds as long 
as Ed/Eo is sufficiently larger than 1, i.e., for 



E 1 < 2.6T (m x /m e ) 2 



(19) 



where Tq = 2.3 x 10~ 4 eV is the current CMB tempera- 
ture. For higher Ej, the integral drops sharply as shown 
in fig. 5. To get a contribution at the lowest Chandra 
energies E 1 = 0.8 keV, for example, would thus require 



m x > 0.6 GeV. 



At energies below that of the CMB photons, Ey / ' z d = 

— 1/2 

2.7 To, (18) scales like E 1 which translates into a sub- 

3/2 

dominant contribution of the form E 1 to the temper- 
ature T A of the excess radio background, where 



T\ 



(20) 



We can put an upper bound on the initial \ abundance 
such that this extra contribution to the radio spectrum 

3/2 

does not exceed the observed E~ i contribution at the 
highest energies where it is reliably observed. Taking this 
to be 10£fl = 4.1 x 10~ 5 eV, corresponding to frequency 
v = 10 GHz, and the corresponding temperature T A = 
0.0032K = 2.8 x 10~ 7 eV (see fig. 1), we obtain the bound 
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on the relative initial abundance of \ to photons, Y x = 
n x /n 7 , before the decays. From (19) one finds that for 
m x greater than a few times m e , the Compton contri- 
bution to the spectrum extends well above the highest 
frequency ~ 90 GeV measured by ARCADE-2 before 
cutting off. For larger values of m x such that the Comp- 
ton contribution extends into the x-rays, we can derive a 
much stronger constraint. 

To compare (18) with observations such as those of 
Chandra [19], we compute 
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85y x (z d /1100)- 1 erg 
cm 2 s deg 

(22) 

in an energy window over which surface brightness has 
been measured (multiplying by 3 x 10~ 4 to convert from 



where the minimum electron energy is E = \E 1 [1 + (1 + str to deg 2 ). The value 85 Y x is for the 0.65 - 1 keV 
ml/EyEy') 1 / 2 ] from the kinematics of Compton scatter- 
ing [32]. To evaluate dn^/dE^ at the present time, we 
must correct for redshift by replacing E 7 — > ZdE~ l and 
E~> — > ZdE^i (where the redshift corresponding to td is 



energy interval, in which the measured value is 10~ 12 
erg cm -2 s _1 deg -2 . We thereby derive a limit on the 
abundance of decaying \ particles relative to photons of 
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Zd — 1) in f c and in Eq, and dividing by an overall 



Y < 10~ 14 



(23) 



which is much more stringent than (21). However we em- 
phasize that (23) applies only for m x > 0.6 GeV; other- 
wise Compton emission in the Ej ~ keV region is kine- 
matically blocked, since then Eo > Ed and the Comp- 
ton spectrum cuts off as shown in fig. 5. In this case, 
a stronger constraint might ostensibly come from lower- 
energy x-ray observations, such as from XMM-Newton 
[20]. But in section 1 we noted that photons of energy 
lower than 1 keV are effectively attenuated through ab- 
sorption by hydrogen. Hence the m x < 600 MeV restric- 
tion is sufficient to generally circumvent the constraints 
from x-ray observations. 




-1.5 -1 -0.5 0.5 



log, n v (GHz) 



6. SYNCHROTRON EMISSION 

To derive the spectrum of synchrotron emission from 
the injected electrons, we follow a similar procedure as for 
Compton in the previous section. The rate of production 
of synchrotron photons from a single electron of energy 
E e is given by [32] 



1 d£L 



dtdE 1 {Ait) 2 v dtdv 

„2^2 

F(x) 



VSirE 2 



(24) 



where x = 2E 1 m 3 e / \3eBE 2 ) and F(x) is the integral over 
harmonics of the gyrofrequency of the electron: 



F(x) 



^3/2 (£R 
47T 2 



3*2*r(|)' 



(25) 



Although F(x) can be expressed exactly in terms of hy- 
pergeometric functions, we find the more tractable ap- 
proximation in the second line, which is just the known 
small- x asymptotic value times e~ x . This is accurate for 
i<1 while overestimating the true value by ~ 20% near 
x = 2.5 where F is maximized, which is sufficient for our 
purposes. 

The desired quantity dn e /dE e is obtained by multiply- 
ing (24) by the electron density n e and integrating with 
respect to time. As in the previous section, the relevant 
time dependence is that of the rapidly decreasing elec- 
tron energy (16). Again changing variables t — > u = x 2 
and taking E 1 — > z^E^ (note that also B — > z 2 d B in or- 
der to reexprcss B at t c i in terms of its comoving value) , 
we find 



FIG. 6: Solid (magenta): the shape of the predicted contri- 
bution to the radio background from synchrotron emission, 
when the constraint (27) is saturated. 



The latter approaches 1.13 as C B — > but falls exponen- 
tially with large C B . To insure that T A ~ E~ 2 - 5 , it is 
necessary to keep C B < 1 over the relevant range of E n . 
This puts a constraint on the injected energy, here from 



m x > 6.7 GeV x 



z d 
1100 



-1/2 



B 



50 nG 



-1/2 



(27) 



The distortion away from the pure —2.5 power law for 
the case where C B = 1, hence (27) is saturated, is shown 
in fig. 6. 

Using T A = (ir 2 / E~f) dn^/dE-y and normalizing to the 
measured value at 45 MHz, we find that the required 
relative abundance is 



Y x = 6x 10" 



z d \ 3 / 2 



1100 



B 



50 nG 



-3/2 



(28) 



which is inconsistent with the Chandra constraint (23) 
unless B exceeds the CMB bound. To make them consis- 
tent, we would need B > 0.8 fiG(z,i/6) 1 ^ 3 - Here we have 
expressed the result with lower redshifts in mind, since 
this reduces the size of the required field somewhat. This 
assumes that m x > 600 MeV so that (23) applies. To 
evade the x-ray constraint by having m x < 600 MeV, (27) 
would require B to be even larger, B > 6.2 /iG(1100/zd)- 
However we shall see that CMB constraints on charged 
particle injection require yet larger values of B than ei- 
ther of these bounds. 



dn 



7,0 



0.048 Y x B 3 I 2 m 



3/2 



dE^ 



^fez^J 2 T Q E'-y 



3/2 



du u 2 ^ 3 e 



u (26) 



c B 



where C B = 2E 7 m 3 e / (3eBz d E 2 ). We are interested in 
parameter values such that the upper limit of the integral 
in (26) is 3> 1; in that case the integral can be approx- 
imated by i the incomplete gamma function r(|,C B ). 



7. DARK MATTER ANNIHILATIONS AND 
NONSUDDEN DECAYS 

In the previous sections, for simplicity we considered 
the case in which electrons are injected suddenly at a 
specific time td (the sudden decay approximation). It is 
straightforward to generalize this to a continuous injec- 



8 



tion of energy; one replaces Y x by 



dt 



(29) 



(taking into account the time-dependence of Zd in the 
expressions for dn 7 /dE 7 ) where dY x /dt is the relative 
abundance per unit time of x particles that decay into 
electrons. The same logic can be applied to annihilations 
XX —> e + e~. The expressions for dY x /dt for decays or 
annihilations are, respectively, 




Y xfi r 1 e */ r , decay 



Y x n 7 (av) , annihilation 



(30) 



Let us consider the synchrotron photons coming from 
annihilation XX ~^ e + e~. We can convert the integral 



over td into one over Zd using td = toz 



-3/2 



which is not 



accurate at late times when dark energy dominates, but 
is adequate during matter domination. Carrying out the 
integral over Zd, where for a first estimate we approximate 
T(5/6, C B ) as a step function that goes to zero at C B = 1, 
we find that dn^^/dE^ in (26) is modified by replacing 



Y z~ 3/2 -> 



^Y xfl t n lfi (av) In 



end 



(31) 



where 2 C nd is the value of Zd such that C B = 1 and 
z dcc = 1100. Here F Xj0 lS the dark matter abundance 
after it has frozen out, which we can rewrite as Y x q = 
n x,o/ n -y,o = 2.7Tof2 x /(m x r2 7 ), and (av) is the cross sec- 
tion to produce e + e~. It is therefore no more than 
(fl DM /n x )(<jv)o, where (av)o = 3 x 10~ 26 cm 3 /s is the 
standard thermal cross section for obtaining the observed 
S] D1I = 0.22. 5 Here we only assume that il x < 0.22 since 
X need not be the only species of dark matter. Now 
equating the modified (26) to the observed radio spec- 
trum, we find 



m x = 0.6 MeV x 



1/2 



B 



50 nG 



3/4 



(32) 



where we define r to be the ratio of (av) for XX ~~ ^ e + e~ 
to the total annihilation cross section. For simplicity 
we assume that any additional annihilation channels do 
not result in charged particles, since otherwise these will 
quickly decay into electrons and behave similarly to the 
electrons produced as primary annihilation products. In 
this estimate we took ln(2dec/^end) ~ 1, which is only 
consistent if z cn( i J$ Zdcc', this puts a second constraint 
on the mass, 



-1/2 



(33) 



5 If there is some other annihilation channel determining the x 
relic density, then {av) for \ — > e + e — can be smaller than this 
bound. 



These constraints begin to overlap (leaving no allowed 
parameter space) if the field is too small, leading to the 
lower bound 



B > 90 nG x 



2/5 



(34) 



In the case where this bound is saturated, the mass is 
m x = 160(rfi x /fi £)A/ ) 1 / 5 MeV. We will see that this is 
strongly ruled out by CMB constraints on charged par- 
ticle injection. Thus dark matter annihilations in the 
early universe cannot explain the radio excess even with 
late-time magnetogenesis. 

For decays, we can carry out the time-averaging to 
find the more accurate dependence of the normalization 
of the spectra upon Zd- This results in a correction factor 
of order unity: 



1 



to/r 



(35) 



where idec = 3.8 x 10 5 y is the age of the universe at 
decoupling, and to = 13.8 Gyr is the present age. For 
tdec C t < (o, the correction factor is approximately 
T(l + 2p/3) = 0.9, 1 for p = 1, 3/2, corresponding to the 
cases of Compton and synchrotron emission, respectively. 
Thus the sudden approximation works quite well for the 
case of decays when r is in the range of interest. 



8. CMB CONSTRAINTS ON CHARGED 
PARTICLE INJECTION 

There are stringent bounds on the injection of ionizing 
energy after decoupling due to its distortion of the CMB 
anisotropies, either through particle decays [17, 33] or 
annihilations [18, 34]. These depend mainly upon the to- 
tal amount of kinetic energy of charged particles that is 
injected, and so are closely correlated with the require- 
ments for producing the observed radio background. We 
first consider the case of decaying particles. Ref. [33] 
obtains 



- <2x 10- 25 s _1 F D 



(36) 



for t > 10 17 s, and a somewhat less stringent limit 
at shorter lifetimes; however the above approximation 
will be adequate for our purposes since the tension of 
this constraint with what is needed for the radio ex- 
cess is minimized at the higher values of r, despite the 
modest weakening of the constraint at low r. Here 
Y DM = (Q DM /Cl^i)(Ej> /m x ) is the standard abundance 
that x would have if it was a stable particle constituting 
the full WMAP DM abundance. Substituting z d for t 
we then find 



Y x < 2 x 10 



-IT 



3/2 



1 GcV 



(37) 
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We have taken Zd — 6 to be the latest decay epoch 
consistent with generating the signal before significant 
clustering of the DM would occur and induce un- 
wanted fluctuations. If we take the synchrotron mass 
bound (27) to be saturated, then (37) becomes Y x < 
10- 18 (6/z d ) 1 / 2 (B/l/iG) 1 / 2 . Comparing this to the re- 
quired value of Y x (28), we get the lower bound on B 



5>5/*Gx 



6 ) 



5/4 



(38) 



For this value of B, the minimum \ mass from (27) 
is m x — 9 x (zd/6)~ 9 / 8 GeV, and the abundance is 
2 x Kr 18 (z d /6)- 3/8 . This is compatible with 



Yx = 



the x-ray constraint (23), which can be written as Y x < 
5 x 10 _17 (Z(j/6). A summary of the bounds on electron 
injection from \ decay is presented in Figure 7. 

For XX e + e~ annihilations, we adapt the result of 
ref. [181 to obtain 




W^H^r'Wo (l^Gev) 



(39) 



where / is an efficiency factor that is close to 1 for anni- 
hilation into electrons, and r was defined below eq. (32). 
Since the rate of injection of charged particles depends 
upon n 2 , we have appended the factor (fl DM /fl x ) 2 for 
the case f2 x < fl DAf to account for the weakening of the 
bound derived by [18], where it was assumed that fl x = 
fl D M- Consider the case where the B field is just barely 
large enough to give the radio excess, so that (34) is satu- 
rated. Then using the fact that (cr v) tot / (&v)o = ^dm/^x 

and the value of m x determined below eq. (34), we find of CS photons, the combination e 4 n 7 Ey 2 coming from 
that the CMB constraint can be written as 



FIG. 7: Summary of constraints on models with electron in- 
jection at redshift Zd from x ~> e + e~ decays, as a function of 
m x . Allowed regions are in directions of arrows. "Isotropy" 
(blue) denotes the limit Zd > 5 due to the isotropy of the sig- 
nal. "Relativistic injection" (magenta) denotes requirement 
limit m x > 2 MeV in order to obtain relativistic electrons. 
The "eq. (27, 28 & 37)" line (brown) is the coincidence of 
upper and lower limits on m x in the case of minimum B field 
given by eq. (38). If B is allowed to be 10 times larger than 
this minimum value, the allowed region is bounded by the 
"eq. (27)" (red) and "eq. (28 & 37)" (black) curves, coming 
from the shape of the synchrotron spectrum and the CMB 
limit on charged particle injection, respectively. 



e 4 n 1 /(E 1 'EQ) is replaced by 



"DA 



< 6 x 10 



-3 



(40) 



A 



2/ x2 "^1/2 , ( ^4 — 1/2 

g (ee) n^E^ 



(ee) n-yE. 



i' 



(41) 



which with eq. (34) then implies B > 3 mG (milligauss) 
in order for annihilations to give the radio excess. This 
is such a strong field that it would overdose the universe. 



9. MILLICHARGED MODEL 

If we insist that B was generated prior to decoupling 
and thus must obey the CMB constraints, additional new 
physics is required. In this section we present a somewhat 
more complicated model, in which the injected particles, 
instead of being electrons, are exotic particles that we 
denote by e, carrying a fractional electric charge ee. In 
adddition, we assume that e has a larger charge g under 
a hidden, unbroken U(l)' gauge group that comes with 
its own "dark photons," (denoted by 7) having a cosmo- 
logical background analogous to the CMB. This scenario 
has been considered previously in the millicharged atomic 
dark matter model of ref. [35, 36]. 

It is straightforward to generalize the Compton result 
(18) to the present model. For the rate of production 



On the other hand in the energy loss rate that eventu- 
ally appears in the denominator, the combination e 4 p 7 is 
replaced by 

A' = g 4 p^ + g 2 (ee) 2 (p^ + p 7 ) + (ee) 4 p 7 (42) 

Thus the CS spectrum will be proportional to A/A' . In 
addition the electron mass is replaced by ms- Assuming 
that the first term in each of (41) and (42) dominates, 
we find that the x-ray constraint (23) generalizes to 



Y x < 5 x 10 




-17 ( z d 



(43) 



In contrast, the CMB constraint depends mainly upon 
the fraction of ionizing electromagnetic energy per de- 
cay. In the present model, this fraction is given by 
(A! — g 4 p^)/A' . The bound (37) is accordingly weakened 
to become 



Y x < 2 x 10 



-17 



8/2 



1 GcV \ 
m x J 



9\ 2 (ft 

[n 



(44) 
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where we have assumed that the leading terms in numer- 
ator and denominator of the reduction factor are dom- 
inant, and also that Tq <C Tq. These two assumptions 
would be consistent if for example Tq ~ O.ITq , the max- 
imum allowed by constraints on the Hubble rate during 
big bang nucleosynthesis. We need to choose the more 
restrictive of the two bounds (43) and (44) in what fol- 
lows. 

The millicharged model can also admit a primordial 
dark magnetic field B whose magnitude is less con- 
strained than that of the visible B field, allowing for 
the possibility of early B magnetogencsis that is uncon- 
strained by the CMB. The modified synchrotron spec- 
trum can be deduced from the derivation of section 6 by 
making the following changes. (1) e — > ee for the rate 
of emission in (24); (2) eB —> gB in the definition of x, 
since the cyclotron frequency is now determined by this 
combination; (3) e — > g in the Thomson cross section and 
p-y — > py coming from the energy loss rate, which is dom- 
inated by 7 Compton emission; (4) obviously m e — > 
everywhere as well. In this way, we find that the required 
abundance (28) generalizes to 

F x , sync = 2.7 x HT 17 x (45) 

which is valid as long as g A Tq > {gee) 2 Tq. Similarly, the 
generalization of the mass constraint (27) is 

^^^(t)" 4 ^)"' (;)'©' (46 > 

Proceeding as in section 8, by inserting (46) into (44) 
and comparing to (45), we find that all the dependence 
upon new parameters cancels out and yields (38) again, 
but now as a constraint on B instead of B. The advan- 
tage is that there are relatively few direct constraints on 
the dark B field, and it could therefore have been gen- 
erated prior to decoupling. For example taking Zd = 6, 
m x = 5 x (m s /m e ) 3 / 2 GeV, B = 5/iG, e = 10~ 3 , g = 
\/T0e, we find from (45) that Y x = 4 x 10~ 12 (ra e /mg) 3 / 2 , 
which marginally satisfies the CMB constraint (44) for 
any value of mg. We can then take mg > 100 MeV 
(hence m x > 14 TeV) and be consistent with constraints 
on millicharged particles [37-39]. The x-ray constraint 
(43) is weaker than (44) and provides no additional in- 
formation. For such values of Y x , the mass density of x 
is sufficiently small before it decays, f2 x ~ 0.04 f^, that 
there is no danger of it significantly affecting the expan- 
sion history of the universe. 

10. CONCLUSIONS 

The diffuse excess radio background has been observed 
by five different groups and therefore appears to be a real 



phenomenon, and the observation of ref. [15] concerning 
its homogeneity seems to constitute a strong argument 
that its origin must be cosmological in nature, in the 
absence of ubiquitous astrophysical sources of size > 1 
Mpc. 

We have shown that injection of relativistic electrons 
after decoupling, through late decays of a metastable par- 
ticle x, can give rise to a radio excess with the correct 
power law fall-off through synchrotron emission. How- 
ever we find that injecting enough charge to get the right 
normalization of the excess background should also have 
distorted the CMB fluctuations relative to their observed 
properties, due to the deposition of electromagnetic en- 
ergy into the intergalactic medium. This contradiction 
can be circumvented if the primordial B field is larger 
than expected on the basis of other CMB constraints, 
B > 5 pG. The size of the required field is smallest when 
the decays are as late as possible (but still before the 
onset of structure formation, to avoid large fluctuations 
in the radio background). Such large fields would have 
to be generated by some unknown cosmological mecha- 
nism after decoupling in order to avoid the CMB i?-field 
constraint, B < several nG. 

Based on our results, the simplest and most conserva- 
tive cosmological scenario for producing the excess radio 
background is through the synchroton radiation emitted 
by electrons from decays \ ~> e + e _ at rcdshift z — 5, in 
a magnetic field of strength B ~ 5/iG, with m x > 9 GeV 
and an initial abundance Y x ~ 10 -18 relative to pho- 
tons. Earlier decay epochs would require larger B fields, 
given by (38). This requires magnetogenesis at low red- 
shifts z < 1100 to circumvent CMB constraints on the 
primordial magnetic field. We find that the annihilation 
scenario, XX ~^ e + e~, is too inefficient and would require 
unacccptably large B fields to produce the observed radio 
excess. 

One might wonder whether decays into other known 
charged particles than the electron could give more 
promising results, since the CMB bound is weakened in 
that case. The weakening comes about because some of 
the initial energy eventually goes into neutrinos, which 
do not ionize the intergalactic medium. We do not expect 
any improvement in this way however, since the same loss 
of efficiency in the production of radio and other electro- 
magnetic backgrounds should occur, and the amount of 
radiation produced by the unstable charged particles be- 
fore they turn into electrons must be very small due to 
their short lifetimes, compared to the time scale for en- 
ergy loss through CS. 

To demonstrate an alternative example that does not 
require late-time magnetogenesis, we presented a model 
of small fractionally charged (~ 10 _3 e) decay prod- 
ucts with mass > 100 MeV, produced from the decays 
of heavy > 14 TeV parent particles. If the fraction- 
ally charged particles undergo cyclotron motion due to 
their interaction with a dark background magnetic field 
B > 5/iG, the resulting synchrotron emission (in normal 
photons) can explain the radio anomaly without conflict- 
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ing with CMB observations. The fractionally charged 
particles are stable relics, but with an abundance less 
than 10 -5 that of baryons, making them innocuous at 
present times. This example might be regarded as an 
existence proof; we hope it will motivate others to search 
for alternative possibilities. 6 
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The integral over xq can be done by expressing 



Appendix A: Stochastic B field from \ 



decay 



In this appendix we estimate the stochastic magnetic 
field spectrum that is generated from the decays of a 
heavy particle \ va ^° e + e~. Consider a decay at time t = 
and position r = 0, producing an electron of velocity 
v. The vector potential due to the e + e~ pair is 



A 



'(\f- vt\~ 



vt 



(Al) 



We can compute the desired correlation function 
(Bi(x)Bj(y)) from (Ai(x)Aj(y)) by taking the appropri- 
ate curls. The correlators can be found by averaging over 
a random ensemble of vector potentials due to currents 
emanating from arbitrary points in arbitrary directions 
(but for simplicity we take them all to arise at the same 
time). We take the contributions from different currents 
to be uncorrelated so that only diagonal terms in the 
product of the two A's contribute. The correlator thus 
becomes a sum over all particles, averaging over the po- 
sition of the decaying parent and the direction of the 
products. Let N be the total number of decaying parti- 
cles in a volume V, and define x± — x — (xq ± vt) where 
xq is the initial position of a fiducial pair; similarly for 
y±. Then 

(Mx)A 3 (y))=e 2 N(v lVj ( fl^l" 1 - (A2) 

x (ly+r-ly-r 1 ))) 

The averaging means 

1 



(X) 



4ttV 



(Txq I dQ v X 



(A3) 



1 



ds 



Jo y/^± 



(A4) 



The integrals over s± can also be done exactly, leading 
to 



(Ai(x)Aj(y)) = e n e j dCl v ViVj (\x-y\- \x+ -y~\) 

(A5) 

where n e is the electron density. The remaining integral 
over angles can also be done exactly, but at large times 
the leading contribution to the gradients will be from the 
first term, which is given by 



4:TTe 2 v 2 n e Si 



p - y\ 



(A6) 



The B field correlator thus goes like l/\x — y\, corre- 
sponding to a spectral index of n B = —2 in the notation 
of section 4. 

The smallest scale that can be relevant for synchrotron 
emission is the Larmor radius r = m e c/eB. Eliminating 
r from this equation and (11) at r = \x — y\, we find the 
field strength B = 47re 3 n e /m e (evaluated at the decay 
rcdshift z d < 1100; this value of B must be divided by z\ 
to convert it to the comoving value for comparison with 
the CMB constraint). Even if the electrons are injected 
with unit abundance relative to photons, this is orders of 
magnitude smaller than the CMB limit and thus irrele- 
vant. 



Appendix B: Brehmsstrahlung 

The final mechanism for photon production from 
charged particles is by Brcmsstrahlung with the ambient 
neutral (before reionization) or charged (after) hydrogen 
atoms. We obtain the photon spectrum in the same way 
as before, by assuming a sudden decay, and replacing 
the integral over time with an integral over the electron's 
energy, which is lost due to Compton scattering. The 
bremsstrahlung spectrum of a relativistic electron collid- 
ing with a hydrogen atom is [32]: 



dA 7 
dEjdt 



[{E 2 e + (E e 



£ 7 ) 2 )0i 



-E e (E e - E y )(f> 2 



(Bl) 



where tq is the classical electron radius and E e — E-y is 
approximately the electron energy after scattering. The 
functions <f>i parametrize the atomic shielding. When 
E 1 <C E e , there are two possibilities: 1) scattering off 
unshielded charges (e.g. in an ionized universe): 



Recently ref. [40] appeared, which also proposes a cosmological 
dark matter mechanism. 



4 In 



E e (E e — En/ 

E 1 m e 



(B2) 



12 



or 2) neutral hydrogen, in which case the shielding goes 
to: 



91 — <P2 



45. 



(B3) 



Then, with E e ^> E 1 we get the result for unshielded 
charges: 



dn 7 2n x nffarQm x 4/6 4 l f m t 



d£\ 



3C e E 7 3 \ m e ' m e ^ _E 7 



8 , / mt 
In 1 



and for shielded charges 

d?i 7 2n x n#ar§m x 4 



di? 7 3C e m e Ey 3 



(B4) 



(B5) 



We have taken the integration limits from m x /2 to m e 
(strictly speaking the lower limit should be m e + _E 7 , 



but the latter correction is negligible at small E 7 ). Note 
that in both cases the spectrum goes as E~ x , which does 
not match the desired shape to reproduce the radio ex- 
cess. We finally estimate the contribution to other limits. 
Rewriting the latter in the language of Section 5, and 
evaluating the spectrum today: 



3an x fiiri e £Ib,o 



dEy AnrripZdE^ Q 



(B6) 



Comparing this with the CS, this is a very small contri- 
bution, even in the keV region: 



-fere 



4 x 10 



-10 



E~, 



1/2 



(B7) 



If injection of electrons instead occurs after reionization, 
the bremsstrahlung contribution is at most larger by a 
factor of 2(3 + 2 ln(m e / 'z d E^)) ~ 80 for z d ~ 10 and v ~ 
GHz, but is further suppressed at high energies due to 
the logarithmic terms. 
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